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Abstract: About twenty years ago, Choptuik studied numerically the gravi- 
O ■ tational collapse (Einstein field equations) of a massless scalar field in spheri- 

cal symmetry, and found strong evidence for a universal, self-similar solution 
at the threshold of black hole formation. We prove rigorously the existence of 
a real analytic solution, that we interpret as the solution observed by Chop- 
tuik. Our construction covers an open neighborhood of the past light cone of 
the singularity. The proof is computer assisted. Starting from an explicit ap- 
proximate solution, we show that nearby there is a true solution. The source 
code and a high precision data file (about 80 significant decimal digits, with 
rigorous error bounds) are included. We do not study perturbations. 
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Figure 1: Scalar field <fi for Choptuik's spacetime. 



x Part of this work was carried out while the authors were visiting NYU Abu Dhabi, in 
fall 2011. 



The Einstein field equations for a four-dimensional metric g with signature 
(— , +, +, +), coupled to a massless scalar field 0, are Ricci^ = 2 d0 ® d0 and 
D g 4> = 0. If one imposes spherical symmetry on g and 0, the problem reduces 
to a two-dimensional problem. The self-similar solution (g, 0) that Choptuik 
[C] observed numerically, and whose existence we prove in this paper, has 
the nontrivial scalar field in Figure 1. 

Figure 1 was generated from an approximation to the true solution, in- 
distinguishable from the true solution at image resolution. Regions with 
< in blue, and with > in red, alternate. The contour lines are at 
= 0, ±|, ±1, ±|, ±2. The self-similarity is captured by a diffeomorphism 

© : Choptuik spacetime — > Choptuik spacetime (1) 

that maps every blue (red) region to the next red (blue) region above it, 
o = —0, and that rescales all lengths by a constant e~ K , equivalently 
Q*g = e~ 2K g. Here K w 1.72 and we prove 

IK-17227262011139106749857559727881918622210 

(2) 
3458805781088634788403570540271455648855- KT 80 1 < 10" 80 v ; 

consistent with the numerical result 3.445452402(3) for 2K in [MG2]. The 
very good review [GM] contains references to other numerical results. 

future singular point (u-,u+) = (0,0) 

central geodesic / \- /""" cr ^ ca ' J ^^™ cone, w_ = 
u _ u ^^ \ /^fs spacetime is analytic across it 

boundary of domain on 
which we prove existence 
w_ = -u + /81 > 



u± < in the 
shaded region 




u_ const: an incoming light cone 
u + const: an outgoing light cone 



Figure 2: The (u_,u + ) plane. The shaded region is the domain of Figure 1. 
At the future singular point, the scalar curvature blows up. Note that the 
dashed line is actually closer to «_ = than shown here. 

We have to specify the coordinates used to generate Figure 1. Every point 
in Figure 1 is a two-sphere. Lines at ±45 degrees are level sets of an ingoing 
or outgoing spherically symmetric null coordinate w_ or u + , respectively. 



The key property that distinguishes the coordinates u± is that is exactly 
a linear reseating of Figures 1 and 2 about the singular point: There is a 
constant \i > such that u± o O = e~ 2nfl u±. This property, and Figure 2, 
determine u± up to a multiplicative constant, see Section 1.4. The constant 
fi « 0.168 is itself a geometric invariant of the Choptuik spacetime. We prove 



I//-1683070789634499695101349790428574207210 



0199080892966476395293134873313662587505- 1(T 80 | < 10" 80 



(3) 



The critical light cone U- = plays a crucial role. For example, by 
causality, any initial data set for the Einstein-scalar-field system that coin- 
cides with Choptuik's spacetime on a spacelike hypersurface E as in Figure 3 
also contains the shaded region above it, and is singular. No such conclusion 
can be drawn for E'. For this reason, it is essential to construct (an open 
neighborhood of) the past light cone of the singularity. 





Figure 3: Principle of causality. 



We prove that (g, (ft) is real analytic on the open domain bounded by the 
dashed line u_ = — u + /Sl > (the value 81 has no special significance). 
This includes real analyticity across u_ = 0, a basic feature of Choptuik's 
spacetime. We now state our result in a logically complete way. 

Main result. Let M = {(«_,«+) G R 2 | u + < 0, u + < u_ < -u + /81}. 
There are constants K, fi > that satisfy (2) and (3), and real analytic 
functions 0, (, Q : M — > R ; that satisfy 
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such that on the four- dimensional manifold M x S , one has /i + Qq > 
with £ = ^ + T^~ , and the spherically symmetric Lorentzian metric 

, , ( 2 (dM- g dn + + dn+ g dn_) £ 2 

5 " exp(-20 ^ /i 2( U+ + U _)2 + ^ + Q£2)2 ^w<^ * 2 

where g stan dard s 2 ^ s ^ e metric on the unit sphere S 2 C R 3 , satisfies: 

• Q*g = e~ 2K g, the spacetime is 'discretely self- similar'. 
Here, is extended trivially to a map M x S 2 — >• M x S 2 . 

• Ricci 5 = 2 d0 ® d0 and D g 4> = 0. 

• The scalar curvature R g = 2g~ 1 (d<p,d<j)) is not identically zero. 
It satisfies R g o = e 2K R g and is therefore unbounded. 

• The boundary u + = U- is a removable standard polar coordinate singular- 
ity, because the functions f = <p,(iQ extend real analytically by reflection, 

f(u-,u+) = /(«+,«_) 

• The solution is close to the high-precision data in the file RefAplusB. dat. 
For a quantitative statement, see equations (27), (44), (60). Also see (61). 

We do not study perturbations about the Choptuik spacetime, that are 
relevant for many of the (conjectured) properties of the solution. Our paper 
can be the starting point for a rigorous investigation of this kind. This would 
be interesting, because current numerical results seem to be inconclusive, see 
[MG1] and [CHLP], and Sections 3.7 and 5.3 of [GM]. 

It has been shown for the Einstein-scalar field system in spherical sym- 
metry that the solution disperses for sufficiently small initial data, and forms 
trapped spheres for sufficiently large initial data, see [Chrl], [Chr2], [Chr3]. 

Since we construct the solution through a contraction mapping, it would 
be possible to state not only an existence, but also a local uniqueness result 
of a technical kind. This would allow us to define the Choptuik spacetime as 
that unique solution. We have taken the freedom to refer to the Choptuik 
spacetime, even though we have not formally stated local uniqueness. 

The rest of Section 1 are detailed (but informal) overviews: 



Section 1.1 
Section 1.2 
Section 1.3 
Section 1.4 



setup and analysis 

basic strategy 

role of the computer 

the invariant fj, and the coordinates (r, £) 




Figure 4: Lines of constant area radius e ^|£|/()U + Q^ 1 



1.1 Overview — setup and analysis 

We reduce the spherically symmetric four-dimensional problem to a two- 
dimensional problem, stated in Section 2.2. Contrary to the introduction, 

we do not use u± as coordinates, and we do not use 0, (, Q as unknowns: 

Coordinates: We use (r, £) given by the change of coordinates 

u a = -(1 + of) exp(-/ir) cr = ± 

Note that (r, £) o = (r + 27T, £), and one can use Fourier series. 

Unknowns: A number /i > (see below) and u = (ui,U2, 003,004) given by 
the change of variables (8). Advantage: The field equations become 
first order and quadratically nonlinear, a property that we consistently 
exploit. Each Ui is 2n periodic or antiperiodic in r, all are 4n periodic. 

The role of \i: For the analysis, \i > is not a constant, but an unknown on 
the same footing as the oo{. This 'gain' of one degree of freedom is, in 
many respects, a compensation for the 'loss' of one degree of freedom 
due to translation invariance of the problem in r. This becomes most 
explicit in (32). 

Note that the ubiquitous, characteristic differential operators (14) are quasi- 
linear but not semilinear, because they depend on the unknown /1. It seems 
that this 'minimal' quasilinearity cannot be removed from the problem - 
not without introducing a complication somewhere else. Fortunately, in the 
space of analytic functions that we work in, the quasilinearity plays a minor 
role. This becomes explicit in the boundedness of the operator in (37b). 
We represent the ooi as combined Fourier-Chebyshev series: 

47r-Fourier series in r — > index m G 7L 

Chebyshev series in £, i.e. 2ir- Fourier series in arccos£ — > index n G Z 
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We work in a space in which the coefficients decay as 0((fti) _ ' m '(K2)~' n ') 
the high frequency limit \m\ + \n\ — > oo, with Ki,k 2 > 1. Such functions 
are analytic on the open subset of C x C in Figure 5, the infinite strip 
| Im(r)| < 21ogKi times the open ellipse with radii |(«2 ± 1/k 2 ). 



±= 



x 



re C 




(eC 



Figure 5: Domain of analyticity. Here ±1 are not the foci. 

We use K\ = 65/64 and n 2 = 5/4. Thus, the domain that we obtain is 
smaller than the one shown in Figure 5. When we restrict to real £, we get 
|£| < 41/40. This explains the earlier inequality u_ < — u + /81. 

The obstruction to using 'big' K\ and k 2 is that one has to avoid sin- 
gularities (e.g. poles) that the complex analytic extension of the Choptuik 
spacetime may have. This motivates our use of Chebyshev series instead of 
power series in £, because a disk with radius bigger than 1, big enough to 
contain the interval [—1,1], may already contain singularities. 

We use a low-high frequency decomposition of the identity, s + (l— s) = 1. 
Here s and 1 — s are projection operators to low and high frequencies, re- 
spectively. In particular, s has a finite dimensional image. These projections 
divide the construction into the computer and analytic estimates parts: 



things related to s — >• 
things related to 1 — s — > 



computer 
analytic estimates 



The shaded region in Figure 6 indicates the image of s. Every • in Figure 
6 stands for a finite number of real degrees of freedom. Only non-negative 
indices m, n are shown, because they determine all, by reality constraints. 
We actually use a much bigger region: 250 x 750. 



m 



" n 



Figure 6: Images of s (shaded) and of 1 — s (non-shaded). 



To work in Fourier- Chebyshev frequency space is very useful, because: 



It reveals the elliptic nature of the problem. In fact, the principal part 
of the equations is morally multiplication by 



im/2 + /i(n + 1] 



(i is the imaginary unit) 



when m, n > 0, with an inverse that is 0((m + n)' 1 ) as m + n — > oo. 
We say morally, because the actual operators in (22), while diagonal in 
m, are only upper triangular in the index n. 

• It disposes of the geometric boundaries (central geodesic and especially 
the critical light cone) in a seamless way, and one does not specify 
initial data anywhere. That's good, because the Choptuik spacetime 
is supposed to be 'universal' after all. By contrast, for alternative 
approaches not based in frequency space, the geometric boundaries 
may be unpleasant to deal with. 

The common slogan in general relativity that 'the constraints are satis- 
fied if they are satisfied initially' does not apply here, because there is no 
initial data in the first place. Not surprisingly, though, the usual differential 
identities yield a second elliptic system, this one linear and homogeneous, 
that is used to show that the constraints are satisfied. All objects related to 
the constraint equations are marked by jj. 

Our analysis uses i 1 norms. This has several advantages: 

• Natural estimate for series convolution. 

• Transparent operator norm. For non- weighted £ l , it is the sup of the 
i 1 norms of the 'columns'. For example, in finite, say two dimensions, 



a b 
c d 



sup 



e 1 <-e 1 



oi,c,(ieR 



If the C 1 is weighted, then the sup has to be inversely weighted. 
And particularly for coding (think finite dimensions here): 

• The i 1 norm can be evaluated over the rationals, for rational compo- 
nents and weights. For this reason, we even use ||^||c — | Re z| + | Ini2;|. 

• If an operator is constructed column by column, its i 1 operator norm 
can be evaluated without ever storing the whole operator at once. 



1.2 Overview — basic strategy 

From an abstract point of view, we construct a nontrivial solution x (think 
metric and scalar field) to a nonlinear system C(x) = 0, where 

C(x) = Ax + B(x,x) 

The unknown x lies in an infinite dimensional real vector space (of Fourier- 
Chebyshev coefficients), A is a linear operator, and B is a symmetric bilinear 
operator. That the Einstein equations, with or without scalar field, can be 
written as a system with only quadratic nonlinearities is a basic fact, but 
often unappreciated. 

Even though the vector spaces are infinite dimensional, imagine that there 
are as many equations as unknowns, i.e. that x and C(x) are in vector spaces 
of 'equal dimensions'. We thereby ignore the constraint equations in this 
overview. The computer is used to construct an approximate nontrivial so- 
lution x rc f, i.e. one for which C(x re f) is 'small'. With x = x re f + x corr , the 
system C(x) = becomes the following system for the correction x co „: 

y/i. "T ^-LJ\%rcfi ' JJ-^corr ^y^teS) -'-'V^'corrj -^corrj \~*) 

with the notation B(a, -)b = B(a,b). The term — £?(x corr ,x corr ) can be at- 
tenuated arbitrarily, by using a sufficiently good x Te f, because it is the only 
term that scales quadratically with the distance of x re f from the hypothetical 
solution x. Somewhat more concretely, one will seek x corr in a ball with small 
radius 1Z > 0, centered at the origin, on which — B(x COTT ,x corr ) will have Lip- 
schitz constant ~ 1Z, and the better x re f, the smaller one can choose TZ, the 
smaller the Lipschitz constant. 

Thus, a good x rc f yields an essentially linear problem for x corr , and the 
key task is to show that the linear operator A + 2B(x m { , • ) has an inverse, 
and control it. We first reduce the task to that of inverting a certain finite 
square matrix X, and a robust method to invert a matrix is given by the 
Neumann series. In fact, if one can find an approximate inverse Y, i.e. such 
that the operator norm ||1 — XY\\ < 1, then the Neumann series 

yj^ L xY- ) = Y ^=o( l - XY ) k ( 5 ) 

is an inverse of X. (It is a right-inverse to begin with, which for finite square 
matrices is also a left- inverse.) This method gives one a lot of flexibility in 
the construction of Y. See Section 1.3 for a detailed discussion. 

The notation used in this overview is local. To connect this abstract 
discussion with the actual construction, we note that the role of C(x) = 
is played by SQ + (fi,u) = 0, cf. (24). With tiny modifications, equation (4) 
becomes the key equation (31). For equation (5), see (58). 
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1.3 Overview — role of the computer 

Computer-assisted analytical proofs have a tradition in dynamical systems, 
beginning with Lanford's [L] in 1982. This paper may be the first in general 
relativity. 

We use the computer for the primitive purpose of number crunching - 
neither for symbolic manipulations, nor for proof verification. Our code is 
portable and yields reproducible results. It uses integer arithmetic only. The 
C source code is in the directory sourcecode (arXiv ancillary files). 

The computer is used for the following tasks: 

(la) To construct a reference (an approximation to Choptuik's solution), 
(lb) To estimate how good the reference is. 

(2a) To construct an approximate inverse (of a finite square matrix). 
(2b) To estimate how good the approximate inverse is. 
(3) To estimate the norm of an operator in Section 4.1. 

From a logical point of view, (la) and (2a) are not strictly required: it 
would suffice to provide the outputs of (la) and (2a) in files. For (la) we 
do exactly that, see the file RefAplusB.dat. For (2a), such a file would be 
quite big, and therefore the code is used to construct the approximate inverse 
column by column, never all of it at once. 

We have included a tool that can generate better and better references, if 
used appropriately, and given enough computing resources. See ImproveRef 
in the file choptuik.c. We have used it to construct RefAplusB.dat. This 
tool is not part of the existence proof. 

From a practical point of view, assuming ~ 2010 equipment: 

(la), (lb) take a few dozen CPU hours, and may be memory intense. 
(2a), (2b), (3) take a few thousand CPU hours. 

The time consuming tasks (2a), (2b), (3) are embarrassingly parallel, they 
split into many individual and completely independent tasks. Every individ- 
ual task can also be run on a personal computer. 

The bottleneck is the convolution of 2-dim arrays of integers, which we 
have implemented using integer multiplication from the GMP library, which 
uses FFT. We use the GCC compiler, and OpenMP. We ran the code on the 
Brutus HPC cluster at ETH Zurich. 

We discuss (2a), (2b) in some detail. Given a real matrix X, the task 
is to construct an approximate inverse Y, see Section 1.2. The matrices are 
N x N, where N = 654375 with our choice of parameters. This number 



arises roughly as 250 • 750 • 2 ■ (0.25 + 0.5 + 0.5 + 0.5) = 656250 with: 

250 number of Fourier indices, < m < 250, see s in Section 1.1 

750 number of Chebyshev indices, < n < 750, see s in Section 1.1 

2 for real and imaginary parts 

0.25 for every field restricted by a Z 2 x Z 2 symmetry 

0.5 for every field restricted by a 1>2 symmetry 

Reality constraints for the zeroth Fourier coefficient yield a value N smaller 
than 656250. This counting applies to the rows of X. For the columns of X, 
one of these degrees of freedom is dropped (a gauge condition that removes 
translational symmetry in r), and one is added for \i. See (32). Of course, 
the role of columns and rows is reversed for Y. 

Each entry of the matrix X is a dyadic rational, i.e. an integer divided 
by a power of 2. The matrix X is not stored entry by entry, but arises 
naturally as a sum of compositions of a small set of operators O. Each O 
is implemented through an efficient algorithm, of which only convolution is 
non-trivial, and each O is equivalent to a matrix with dyadic rational entries. 

There is no useful linear order for the rows and columns of X, because 
there are two frequency indices m and n. This makes it hard to visualize X. 
We will simplify our discussion below by pretending, very informally, that 
there is just one frequency index, and that the rows and columns are arranged 
in order of 'increasing frequency'. With this simplification, the matrix X has 
the structure in Figure 7. 




Figure 7: Matrix X, with N' w iV/100. 

The square submatrix O is dense and some of its entries are large. The 
band © is close to the identity, but not very close in the operator norm. The 
entries in © are extremely small in the operator norm. The entries in O are 
identically zero. To summarize, X is close to a band matrix in the operator 
norm. However, the width of the band is in the thousands. 
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In (2a), the construction of Y, only O and ® are taken into account. We 
first apply Gaussian elimination to O, with partial pivoting and rounding to 
dyadic rationals at each step. The output of Gaussian elimination is com- 
bined with @ to a first naive approximate inverse Y naive . Finally, a truncated 
Neumann series is used to generate from Y naivc a better approximate inverse 
Y. (For efficiency, our implementation differs slightly from this description.) 

This approach to (2a) exploits the approximate sparseness of X, and it 
made the calculation feasible with the resources available to us. In principle, 
with enough resources, one could apply Gaussian elimination to all of X at 
once. 

In (2b) we calculate the operator norms ||y|| and ||1 — XF||. This is the 
only place where © has to be used. Only once we have checked \\1—XY\\ < 1 
do we know that the construction in (2a) was actually successful. We can then 
use the Neumann series inequality ||X _1 || < ||y||/(l — ||1 — XY"||). 

Tasks (2a) and (2b) are embarrassingly parallel, because both the con- 
struction of Y and the evaluation of ||y|| and ||1 — XY\\ can be done column 
by column, see the discussion of the i 1 operator norm in Section 1.1. 

Tasks (2a) and (2b) have to be repeated for the Jj system, that is used to 
show that the constraints vanish. Here, N = 280125 « 250-750-2- (0.25+0.5). 

1.4 Overview — the invariant /i and the coordinates (r, £) 

The purpose of this section is to motivate the geometric ansatz that we use to 
construct the Choptuik spacetime. It provides important intuition. However, 
it is not part of the existence theorem. First, we need three facts. 

Let G be the set of germs of C°° functions R — > R at the origin x = 
0, with /(0) = and /'(0) > 0. Germs are convenient for making local 
statements about functions, without having to worry about domains; for a 
definition see Wikipedia. The set G is a group under composition (o). Two 
elements f,g G G are called conjugate if and only if there is a n G G such 
that k o / = g o n. For all f,g,KE G: 

Fact 1. If f,g are conjugate, then f'(0) = g'(0). 

Fact 2. If f'(0) y^ 1, then f is conjugate to the linear function x h- >■ f'(0)x. 

Fact 3. // / is linear with /'(0) ^ 1, and f o k = k o f , then k is linear. 

In words: The derivative at the origin is a conjugation invariant (Fact 1), 
and if it is not equal to one, then it is the only invariant (Fact 2). 

The rest of this Section 1.4 is not rigorous. We blur the distinction 
between a function and its germ, and we use implicitly some very general 
qualitative properties of the Choptuik spacetime. 
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The Choptuik spacetime has a conjugation invariant as in Fact 1. To 
see this, first introduce a set of incoming null coordinates. Namely, let 
NullCoorcL be the set of all u : (Choptuik spacetime) — > K, for which: 

u is a spherically symmetric incoming null coordinate, with du ^ 
u = is the critical light cone of the Choptuik spacetime 
u < in the interior, u > in the exterior 

The group G acts on NullCoord_, in fact 

G x NullCoord_ -> NullCoord_ (/» H- / o <u 

is well defined as a map, and is a group action. This action has an important 
property: for any two u, u' in NullCoord_, there exists a unique /eG with 
/ o u — u' . Henceforth, this property is called regularity. 

Let be the self- similarity diffeomorphism (1) of the Choptuik spacetime. 
If u is in NullCoord_, so is u o -1 . By regularity, there is a map 

NullCoord_ — > G, u \—t f u defined by u o 6 _1 = f u ° u 

We have f K0U = K o f u o k~ x for all kgG. By regularity: 

By Fact 1. /^(0) is an invariant (the same for every u e NullCoord_). 

The constant fi used in this paper is defined by e 2 ""^ = /„(0). We anticipate 
that for the Choptuik spacetime, /i > and /^(0) > 1. Therefore: 

By Fact 2. There exist u e NullCoord_ such that /„ is linear. 

By Fact 3. If /„ and /„/ are linear, then /„ = f u i and u = cu' for a c > 0. 

Fix aw_6 NullCoord_ for which f u _ is linear, u_o9 _1 = e 27rM w_. Let u + 
be the spherically symmetric outgoing null coordinate that satisfies u + = u^ 
along the central geodesic. Then u + o © _1 = e 2nfl u + , because it is true on 
the central geodesic. The 2D coordinate pair (u-,u+) is unique up to joint 
multiplication by a positive constant. 

The 2D coordinate pair (r, £) used in this paper can now be defined by 

u a = -(1 + cr£) exp(-yur) a — ± (6) 

By construction, £ = is the central geodesic, £ = 1 is the critical light cone, 
and (r, ()o0 = (r + 2n, £). The pair (r, £) is unique up to adding a constant 
to r. We compare the two coordinate systems in a table: 





Self-similarity diffeomorphism G 

(with o = —0) 


Null directions 
(2D causal structure) 


(m_,u+) 


(u_,u+) o 6 = e~ 27r ^(M„,M + ) 


ker du a 


(r,0 


(r,Ooe = (r + 27r,0 


ker[crd£-/i(l + (xO dr ] 
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In the first case, depends on /x, but the 2D causal structure doesn't. 
In the second case, doesn't depend on /x, but the 2D causal structure does. 
Either way, the constant /x > appears somewhere. Since it is a geometric 
invariant of the Choptuik spacetime, /x cannot be chosen freely before the 
construction. Rather, it is an unknown and it has to be constructed. 

2 Einstein equations coupled to scalar field 

The equations for a metric g, coupled to a massless scalar field 0, are 

Ricci 9 = k d<p <S> d<ft Einstein field equations 

n g 4> = wave equation 

where Ricci 9 is the Ricci curvature of g, the one form d0 is the differential of 
0, and D g is the Laplace-Beltrami operator for g. All values for the constant 
k > are equivalent, by rescaling 0. This rescaling is necessary to compare 
papers that use different values for k. This paper uses k = 2. 

The four- dimensional problem reduces, in the case of spherical symmetry, 
to a two-dimensional problem: 

Section 2.1 




abstract 4D problem 




From a pedagogical perspective, the direction from 4D to 2D (Section 2.1) 
is easier to understand, because it starts from the basic equations. From a 
logical perspective, only the direction from 2D to 4D (Section 2.2) is required 
for the existence theorem that we prove. 

We use a frame formalism, because it yields first order, quadratically 
nonlinear equations. The 2D formulation has more equations than unknowns. 
The equations are dependent, see the differential identities in Section 2.3. 

2.1 From 4D to 2D 

Let (r, x 1 , x 2 , x 3 ) be Cartesian coordinates on R 4 . Spherical symmetry means 
invariance under the standard action of the matrix group 0(3) on the co- 
ordinates (x l ,x 2 ,x 3 ), with t fixed. It is convenient to start with Cartesian 
coordinates, because (as opposed to polar coordinates) they are regular along 
the line (x 1 , x 2 , x 3 ) = 0, the set of fixed points of the action. Set 

f = v^x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 
13 



Let /i > be a constant, and let Q and ( be real valued functions that 
depend only on r and £. Introduce a frame (e ,e!, e 2 , e 3 ) and a spherically 
symmetric metric g with signature (— , +, +, +): 



d_ 



■\- n^x 



k=\ 



dx k 











fc=l 



, d 



x 



dx k 



with i — 1, 2,3 



(7) 



# x = exp(2() ( - e ® e + e x <g> e : + e 2 !g> e 2 + e 3 ® e £ 

A spherically symmetric spacetime can always be brought into this form, by 
appropriate choice of the two free functions Q and (. We omit a more precise 
statement or a proof, because it is logically unnecessary for this paper. 
We assume /i + t^ 2 Q > 0, which implies that the frame is non-degenerate: 



e A ei A e 2 A e 3 = n(n + £ 2 Qf 



d d 

A 



A 







A 







dr dx 1 dx 2 dx 3 



By construction, the functions u a = —(1 + cr£) exp(— fir) with a — ± 
both solve the Eikonal equation g~ 1 (du a , du a ) = 0. This determines the 
causal structure of the restriction of the metric to the 2-dimensional plane 
( r ;0 l—J * ( r ?^0,0). This causal structure depends only on the constant //, 
not on Q or (. See Figure 8. 



*r 


/ \y\^ i^ 


/ \/Y/ l^ 









1 



s c 



Figure 8: Forward light cones. Exact directions depend on \i. 

However, £ = 1 is a null hypersurface for all /i, because it is a level set of u. 



Introduce the two radial null vector fields 

x k e k 
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D° = ae + J2 



In particular D a {r) = o and D a (£) = /x(l + cr£). Let <fi be another function 
that depends only on r and £, it plays the role of the scalar field coupled to 
the Einstein equations. The functions £, Q, <fi are not themselves useful for 
us. Instead, the basic unknown functions of r and £ £/iat are used throughout 
the paper (together with the unknown constant \x) are 

wi = 2^g + d- (c + io g (/i + e 2 g)) + ^ + (c + io g (/x + eQ)) 

ul = D°{C) -oil 
< = D^cf)) 

with a — ±. A useful property is the invariance of these seven functions 
under a global conformal transformation of the metric, £ — > £ + const. 
The equations (8) can be rewritten as 

g=( Wl -o;J-w+)/(20 (9a) 

Zr(log( Ai + Qa)=^2-^3 (9b) 

D'(C)=wf + (r/i (9c) 

D°{<l>) = ujI (9d) 

This can be used to calculate Ricci 9 (e m ,e n ) — 2e m (0)e„(0) with m,n — 
0, 1, 2, 3 and exp(— 2£)n 9 in terms of just u>i, wf, ^3, wj. This lengthy but 
straightforward calculation is not included. The result is 

x 

Ricci 5 (ej, e,-) - 2e i ((j))e j ((f)) = -ri{^i + ^1" + ^ + ^ - ^ _ ^2tt) 

+ ^- (ftj- + fi+ - fi^ - fi+ + n^j, + n+| + 2Qg + 2fi+) (10a) 



for all z, j = 1, 2, 3 and 

1 

2l 



Ricci 5 (e , e ) - 2e o (0)e o (0) = — ( - Q 2 - Q$ + fi 2fl + ^J - ^3 - ^3") 



Ricci 9 (e , e<) - 2eo(</>)e;(0) = — ( - Q x + ft+ - 0, + fi+) 

exp(-2C)n 9 = — (fi 4 " + fi+) 



(10b) 



Here, by definition: 

nf(//, w) = £(£" + ff/i)(wi) + A*(2wi) 

+ ^(cuiCJi — 2co'2 _cr (jJi — 6u%(jJi + 4wjWi + u; 2 ~ CT u; 2 ~' J (H a ) 



O", ,c^ 



- 2w 2 T o; 2 " <T - 4w> 2 - <T + o;£a;£ + 4w£wf - 4u%u 
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;iib) 



lie) 



fi£(/*,w) = £(£>" + a/i)(cJ 2 " CT ) + ju(wi + wj ff + wf) 

+ ^(wiWi — 2ix>2~Wi — Iuj^Ui + UJ2U2 ~ 6^2" w^" + wjo;^) 
fi£„0i,w) = £(£>* + ^)K) + MM - 2o£) 

+ it(-4u>5uZ + 8^ - 4«) 
fi£(/x, w) = f (£T + a/i)K ff ) + M~Wi) 

+ |£(— u)iU)i + 2dJ2"a;i + 20J2OJ1 — OJ2OJ2 + ^£^2 (Hd) 

— UJ2OJ2 — ^oj^uj^) 
fij(^,w) = Z(D° + o-MK"*) + /i(w 4 - ff + u£) 

+ ie(-4w 4 + w 2 " - 4w 4 -w 2 + ) 



(lie) 



Claim. // Ricci 5 = 2 d.0 <E> d<p and n g <f) = 0, then the ten functions VL? with 
i G {1,2, 2ft, 3, 4} and a — ± a// vanish. 

In fact, [D - , L>+] = M-° _ + D + ) and equations (8) imply that VL<1 - tt* - ^ 
and — Q2 +^2" an d — fi^ +^3" and —Q£+Q~[ all vanish. By Ricci 9 = 2 d0®d0 
and D g (f) = 0, all $Ts vanish. 

Remark. The principal parts of the right hand sides of (11) have the di- 
rections indicated in Figure 8, and their coefficients vanish at £ = 0. This 
degeneration at £ = is standard for polar coordinates. 

2.2 From 2D to 4D 

A periodic 2D problem (called 2Dprob) is defined below, every solution of 
which yields a 4D spacetime as in Section 2.1. The periodicity of 2Dprob 
corresponds to the fundamental self-similarity diffeomorphism G of the 4D 
spacetime. This section is logically independent from Section 2.1. 

Prerequisites (2Dpre). Let £* > 1 be constant. Let 

(r,0 6 3 = (R/4ttZ)x (-£*,£*) (12) 

The set 3 is a cylinder. The unknowns are a constant \i > and four real 
analytic functions Ui,U2,U)3,U4 : 3 — > R that satisfy (a) through (d): 

(a) Under (r, £) \-¥ (r -+- 27r,£) ; they transform as 

(wi, w 2 , w 3 , W 4 ) I-)- (Wi, w 2 , w 3 , -w 4 ) 
TTiai is, i/iree functions are invariant, one changes its sign. 

(b) If 0J4 is invariant under (r, £) H- (r + T, £), i/ien T e 47rZ. 
In particular, a; 4 does no£ vanish identically. 
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(c) The function oo\ satisfies Pou\ = —ooi. See (13). 

(d) The inequality 2/x + £(u>i — oo 2 + Pw 2 ) > holds on 3- 5ee (13). 
Here and below, P is defined by 

(Pf)(r,0 = f(r,-0 (13) 

for every function f . Set 

D* = a£ + H(l + a£)% a = ± (14) 

They are well defined as vector fields on 3- Set u^ = Ui and oof = —Poji for 
i = 2,3,4. Adopt equations (11) as definitions. The Q's are now real valued 
functions on 3- 

Periodic 2D problem (2Dprob). Find fj, and (001,002,003,004) that satisfy 
2Dpre and that satisfy Qf = £7^ = Qi = ^3" = ^4" = identically on 3- 

Claim. Every solution ft, (001,002,003,004) to 2Dprob yields a unique real 
analytic solution (g,4>) to Ricci 9 = 2d0(g)d</> and D g (f) = 0, as in Section 2.1, 
on the open set R x B^(0) C R x R 3 ; that satisfies: 

• The constants /i in Sections 2.1 and 2.2 are identified, and oo\, oof, oof , 
oof in Sections 2.1 and 2.2 are identified as functions ofr,^ for £ > 0. 

• ((0, 0) = and </>(2tt, 0) = -0(0, 0). 

Furthermore, this unique solution satisfies for all (r,x) 6 R, x B^(0): 

Qo&= Q (15a) 

(o9= ( + K with K = 27ifx- J^(1tuj 3 (t,0) (15b) 

0o0 = -0 (15c) 

where : (r, x) *— > (r + 2ii,x) is a global conformal isometry of g. 

To check this, suppose /i, (001,002,003,004) solve 2Dprob. Note that PD° ' f = 
—D~ a Pf for every function /. Recall that P001 = —001 and Poof = —oo^ a 
with i = 2,3,4, by (c). Therefore, equations (11) imply PQ° = —tt~ ff . It 
now follows from 2Dprob that Q~[ = Q2 = ^Vh = ^3 = ^4 = on 3- That 
is, all Q's vanish on 3- With this piece of information, we can introduce 
two real analytic functions Q, on 3, and one real analytic function ( on 
R x (—£*,£*), the universal cover of 3, as follows: 
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• Let Q : 3 — > R be given by (9a). It satisfies (15a) and PQ = Q. 

• Let ( : Rx (— £*,£*) — > Rbe the unique solution to (9c) and C(0, 0) = 0. 
Such a solution exists, because the one form a on 3 given by a(D a ) = 
Co>3 + aji is closed: 

£(da)(Pr, D + ) = £Pr(a(P>+)) - £P>+(a(Pr)) - £a([£r, D + }) 

= t(D- - ri(tut) - t(D + + rifa) = «s -nt = o 

Thus, da vanishes when £ ^ 0, and by continuity at £ = 0. Then 
f = (oQ— ( satisfies d/ = 0, and is equal to a constant K. The equation 
(JK)(t,0) = /i - w 3 (t,0) implies (15b). The difference f = P( - ( 
satisfies d/ = and /(r, 0) = 0, and therefore P£ = ( identically. 

• Let : Rx (— £*, £*) — >• R be the unique solution to (9d) and 0(27r, 0) = 
—0(0, 0). Such a solution exists, because the one form a(D a ) = ul is 
closed by fij — fij = 0. The sum / = o + satisfies d/ = and 
/(0, 0) = 0, which implies (15c). It follows that 0(r + 4vr, £) = 0(r, 0, 
i.e. is a function on the cylinder, : 3 — > R- Finally, P0 = 0. 

Denote by the same symbols Q, C an d the corresponding functions on R x 
-B^(O) C R x R 3 . Define a frame (eo,ei, 62,63) and a metric g by (7). The 
condition \x + £ 2 Q > follows from (d). The frame, the metric and are 
real analytic (here PQ = Q and P( = £ and P0 = are used). While (9a), 
(9c), (9d) hold by definition of Q, (, 0, equation (9b) is now a consequence of 
equation (9a), of the fact that the definitions of D ± in Sections 2.1 and 2.2 
are consistent, and of Q^ — Q 2 — fi^ = 0. Therefore, equations (8) hold. The 
vanishing of all the fi's and equations (10) imply that Ricci 9 — 2 d0 (g) d0 = 
and D g 4> = when £ 7^ 0, and when £ = by continuity of the left hand 
sides. The uniqueness part of the claim is automatic, because every step of 
the construction is forced. 



2.3 The jj system 

Suppose fi, (cui,cu2,u 3 ,u>4) satisfy 2Dpre, but not necessarily 2Dprob. A 
straightforward calculation with [D~ , D + ] = fi(D~+D + ) and — D~ = PD + P 
shows that the £Ts satisfy a homogeneous differential identity. Namely, with 
Vt-i = Q~ for alH G {1, 2, 2JJ, 3, 4} (equation may continue on the next page): 

i(l-P)( J D+ + /i)fi 1 \ f |(1 + PR 



\(D +2 + fi)Q 2i - P(D + + v)VL 2 ) + ^ V^i I" Pfi 2 - 2P0 3 
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£ f f2Pu 2 - Pu 3 +ui — 2u 2 + co 3 \ fix + PQi 



2 \\ Pw 2 + uj\ - uj 2 J 2 

P0J 2 — P^Z + w 2 ~ ^3 \ ^1 — -P^l 



Pw 2 + Wi - u 2 J 2 

—P00 3 — 3cJi + CJ 3 \ f2 2 + Pfi 2 



-Puj 2 — uji — 7(jj 2 + 4:Cj 3 J 2 

-P0J2 — P0O3 — 0O 2 — U> 3 \ ^2 — P&2 



■Pu 2 — Wi + w 2 — 4w 3 y 2 

/ Pcj 3 — cui — u 3 \ 2 (j + Pf2 2t j 
y3Po;2 - Wi + w 2 y 2 

Pw 2 + Pw 3 + u 2 + w 3 \ fi 2 ft — Pfi 



«2j 



3Pw 2 - Wi + u 2 J 2 

2wi\ fi 3 + Pfi 3 f-2Pu 2 - 2oo 2 \ fi 3 - Pfi 3 



4w 2 y 2 V _4w 2 / 2 

2Pw 4 - 2w 4 \ fi 4 + P^4 /2Pw 4 + 2w 4 \ fi 4 - Pfi 4 



-4w 4 ) 2 \ 4w 4 y 2 

Claim. Suppose n and (oo\, 002,003,004) satisfy 2Dpre. // 

(1 + P)0+ = 0+ = fi+ = fi+ = on 3 
t/ien Pfii = — f2i and ('i/us will be referred to as the Jj system): 

, f f /Pa; 2 - Pw 3 + w 2 - w 3 \ q , f P^s -ooi~oo 3 \ Ogj + P^2(t -gx 
2 \ V Pu 2 + U)\ - uj 2 J x \3Pcu2 -001 + 002) ~ 2 { ' 

P0O 2 + PUJ 3 + U>2 + CJ 3 ^\ ^2() — P^2i 
3PU2 — 00 1 + UJ2 

3 Fourier- Chebyshev series 

We reformulate 2Dprob in Section 2.2 as a new problem 2DprobSeries for 
Fourier- Chebyshev series, see page 27, and introduce handy notation for all 
the objects in 2DprobSeries, to prepare for the analysis in Section 4. 

This section is essentially a collection of definitions. Reading through them 
linearly makes sense from a logical point of view, but is not advisable, because 
the motivation for each definition lies in the definitions that come after it. 
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Series. For every v = (v mn ) in the real vector space 

J v mn = v_ m _ n and v mn = v m - n and there 

v — \ {Vmn)m,n£W, C (L , , , , 

l area,/? > with sup minez |e Q l m l +/3 l"l^ mn | < cx). 

we introduce the Fourier in r, Chebyshev in £ series 

ToSeries (w) = ( £) m nGZ v mn exp(±imr + in^) ) (17) 

It is a function of (r, £), or a function of (r, 6*) that is even in 9. 

Properties of (17). Let v E V. 

• There are a, j3 > as above, depending on v, such that (17) converges 
and is analytic jointly in (r, £) on the open, product subset of C x C: 

|Im(r)l<2a and (c^ + (kdw <x 

The ellipse contains the real interval £ G [—1,1]. Cf. Figure 5. 

• The series (17) is real when r and £ are real. 

• In terms of the polynomials Chebyshev n (£) = cos(narccos£), 

ToSeries(u) 

= Emez ^moexp (\imT) + 2E m62 , fl> o%nexp(ijmr) Chebyshev n (£) 



Convention: In this section, ToSeries is thought of as a map from V 
to maps (R/47tZ) x [—1, 1] — y R,. Advantage: all functions produced by 
ToSeries have the same domain. By analytic continuation, the identities 
that we state will automatically hold on appropriate domains inCxC. 



Fields and multifields. We refer to elements of V as fields, and to tuples 
of fields as multifields. For fields, we define the real vector spaces: 

VWrwoPiPeriodic = {v E V \ v mn = for all odd m} 
VxauTwoPiAntiPcriodic = {v E V \ v mn = for all even m} 
Vcaugcd = {vEV\ Re(v 10 ) = 0} 
Vxiodd = {v E V I v mn = for all even n} 
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For multifields: 

W = V TP x V T p x V T p x V T p x V T ap C V 5 

V = (V TP n Vxiodd) x V TP x V TP x V T ap C V 4 

Vcaugcd = (V TP n VxiOdd) X V TP X V TP X (VpAP n Vcauged) C V 4 

V«=(V T pnVxiOdd)xV T p cv 2 

With the local notation VtP = VpauTwoPiPeriodic, VpAP = VrauTwoPiAntiPeriodic- 



Every field v G V is determined by its m, n > elements. Precisely, the map 
f _ v o n G R and there are a, (3 > 

V — V \ {Vmnjm n>0 (— ^ _i_a 

[ with sup mn > \e am+pn v mn \ < oo 

is bijective, an isomorphism. This useful fact will often be used implicitly. 



Auxiliary linear operator U. Let / = {f mn )m,n>o C C be a sequence with 
fo n G R, and for which \f mn \ has a polynomial bound in m and n. For k > 
define U [kj] : V ->■ V by 



Vm, n > : (U[A;, f\v) mn = f mn v m , n +k 

It is a well defined V — >• V map. It acts diagonally on the index m, and there- 
fore leaves the subspaces V Tau TwoPiPcriodic and V Tau TwoPiAntiPcriodic invariant. 

As a special case, introduce the operator Uff x : V — )■ V by 

if n = 



Vm, n > : (U.f ^ 



urn 



w m , n _i if n > 



Auxiliary linear operator V. Let / be as above, with \f mn \ < 1 for all 

m, n. For k > 0, define the linear map V[/c, /] : V — > V by 

Vm,n > : (V[fc, /]v) mn = EZoi-^iKf^mn 

= Ei=o( U [ fcL '/i] ?, )™n 
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with (f L )mn = (-^) L U.£ = d fm,n+u- Note that 

V[kJ}(l + XJ[kJ]) = (l + \J[k,f])V[kJ} = l 
In particular, 1 + U[k, f] is invertible if k > and \f mn \ < 1 for all m, n. 



To see that V[k,f] : V — > V is well defined, note that K/zOmnl < 1 for all L,m,n. 

:> m,n>0 I 



Therefore, if a, /3 > are such that sup m „ >0 |e am+ ^"u mn | < oo, then 



sup m ,„>o|e Qm+/3 "(V[fc,/] W ) m „| < (1 -e-^)- 1 S u Pm ^ \e am +^v mn \ < oo 
In other words, if a, j3 > work for u, then the same a, /3 also work for V[k, f]v. 

Notation for U, V. Three examples illustrate the notation. Here i = \/—l. 

V[2,-l] =V 2j _! =V[2,/]with/ ron = -l 

U[0, <5 0re ] = U ,5o„ = U[0, /] with f mn = 5 0n (Kronecker delta) 

U[l, -im/2 + fm] = U x _ im/2+Mri = U[l, /] with f mn = -im/2 + \xn 

Basic linear operators. Introduce the V — > V linear operators 

P — U 0] (_l)n 

TauDerivative = U ,j m /2 

E=i(U5 il + U 1 ,i + U llJbn ) 
OnePlusXiTimesXiDerivative = (2Vi ) _i— l) Uo, n 

S-^eg = 2 U 1 , 1 V 2jl 

They satisfy (here P S ection 2.2 is (13)) 

ToSeries(Pu) = P Section 2.2 ToSeries(w) (18a) 

ToSeries (TauDerivative v) = §p ToSeries(t> ) (18b) 

ToSeries(H v) = £ ToSeries(u) (18c) 

ToSeries(OnePlusXiTimesXiDerivative v) = (1 + £)^ ToSeries(t;) (18d) 

ToSeries(H- 1 - rcg tO = ^(ToSeries^;) 

-ToSeries (^l^o) (18e) 
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The operator S ' rcg is only a left-inverse of 



SS-^ = l-(U 0An -UyjV 2il (19b) 



The identities (18a), (18b) are checked directly, and for (18c) one only has to use £ = 
cos(9 = \{e ie + e~ ie ). To check (18d) and (18e), observe that 

sin6^ToSeries(t>) = ToSeries(±(U!f A - U M - U Mo J U ,„ v) (20a) 

TbSeries(t;)| c _M) = ToSeries((Uo,«5 „ -U 2 ,5 0n ) V 2 ,i«) (20b) 

with £ = cos 6* in (20a). Here sin 6*^ has to be moved inside ToSeries as a whole, because 
the definition of V only allows even functions of 6. Multiply (18d) by (1 — £) and use 
(1 -<f )j| = -sinfl^ and (20a) and (18c). Multiply (18e) by £ and use (18c) and (20b). 
This reduces (18d) and (18e) to the operator identities 

(U?,i + Ui.i + U Mo „ -2)(2 Vi.-i -1) U ,„ = (Iff, 1 - Ui.i - U Mo „ ) U 0> „ 
(U?,i + Ui,i + Ui, 5o J Ui, x V 2 ,i = 1 - (Uo, 5o „ - U 2A JV 2]1 

The next step is to remove all Vs. In the second line, just multiply from the right by the 
invertiblc 1 +U2 1 to remove V 2 1. The first line requires an extra step: add the vanishing 
term 2 Uo,5„ Uo, n to the left hand side, obtain a sufficient identity by removing the three 
occurrences of Uo, n , only then multiply from the right by the invertiblc 1 + Ui,_i = 1 — Ui,i 
to remove Vi,_i- Check the resulting quadratic U identities using Uff 1 Ui,i = 1 — Uo : 5„ 

and Vk,fVe,g = Uk+£,h with h mn = fmn9m,n+k- 



Convolution. Define the R-bilinear convolution * : V x V — > V by 

Vm,n G IL : [V * W) mn = 2^ mi,ni,m 2 ,n 2 eS: v mmi w m2n2 

mi+m2=m, ni+n2=n 

It satisfies ToSeries(w * w) = ToSeries(f) ToSeries(u>). 



The J and K operators. Let \x > 0. Introduce the V — > V linear operators 

/ = OnePlusXiTimesXiDerivative +1 
J^ = TauDerivative +/i /f 
/ K=I(J+P / JP) 

K^ = TauDerivative +/i/K = ^(J M +P J M P) 
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with ToSeries(J M u) = (P^ ection 22 + /i) ToSeries(w). Note that 
1(1 + crP) J ±(1 + <tP)=K 1(1 + aP) a = ± 



(21) 



arc: 



Equivalent definitions of the above operators, and expressions for the inverses 
of J and K, with L M = U , im /2+M™+i) and ^ = U t ' im/2+]i( n +i) 1 ' 

4 = Vi_i(U ,n+l + Ui, n ) 



and 



J/u — Vi : _i ( Uo,j m /2+/x(n+l) + ^l,-im/2+nn ) 

= V W L M (1 + U[1, T -J^g_]) 

j; 1 = v[i, . H ^SSi) ]V(i+ui.-i) 

/ K = V 2 ,_ 1 (U ,„ + i + U 2 , n+1 ) 

K M = V2,-l ( Uo,+j m /2+^(n+l) + U2 j _j m /2+p(n+l) 



K,- 1 



m ::ia;;! ]V(i+vi) 



(22a) 



(22b) 



To check ^t in (22b), note that PV k j P = V k 



and PVkjP 



V fc,(-i)*/) and 



>k,(-l)Xf 

2ft =fi +P/1 P = Vi,_i(U 0l „+i + Ui, n ) + (P Vi,_i P)(P(U ,„+i + Ui,„)P) 
= Vi 1 _i(U ,„+i + Ui,„) + Vi,i(U ,„+i - Ui,„) 
= V 2| _i(l + Ui,i)(U ,„+i + Ui, n ) + V 2 ,_i(l + Ui,_i)(U 0l ^.i - Ui, n ) 

Last step: V2,-i = Vi l<T Vi,-^ for cr = ±, because 1 + U2,-i = (1 + Ui,_ CT )(l + Ui i<T ). 



Linear operators for multifields. Set 



■=:-!, reg 



/id 



V 



/f(l + P) 0\ 

10 

10 

\ 1/ 

P) 0\ 

10 

10 

1/ 



o c 



l,reg 



P) 




-POO 



: W^ V 



W^ V 
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The operators S and S" will be used to split the equations into two parts. The 
part selected by S has 'as many equations as unknowns', the part selected 
by S" are the constraint equations. For fi > introduce 





(^ o 0\ 




J At 


H,= 


-J^P 




J M 




{o o jj 


fS = di&gifiJJJ) 


J M = diag(K M , J M , J M , J M ) 


= TauDerivative +fi fl 



v^w 



V^ V 



V^ V 



and note that J^ 1 = diag(K *, J 1 , J x , J^ 1 ). Note that/T, J M map V — > V by 
[P„K] = [P, K M ] = 0. The definition of V and (19a), (21) imply J M = S E H^. 



The operators r\ and T 2 . Define the linear operator Vi : V — > W by 



/ 2 Vl 

Ul + t> 2 - Pt'2 

T lV = 2Pw 3 - 2Pt; 2 
-vi 
\ Vi~ Pv 4 J 



with v = (vi, . . . , 1)4). As required, r\v G W. Define the symmetric, bilinear 
operator r 2 : V x V — > W by (equation may continue on the next page) 



To(v,w) 



1 
4 



+ 



1 



fSPv 2 - 2Pv 3 + vi- v 2 \ 

Pv 2 +V1 — V2 

* w\ 

-Pv 2 — Vi + v 2 

(PV 2 + 2PU 3 + V! + v 2 \ 

2Pv 2 + 2v 2 

-2Pv 2 + 2Pt; 3 

-Pv 2 - v 2 

\ PV4 + V4 ) 



w 2 + Pw 2 



( ~2v x 

Pv 2 -V\-V 2 

2Pv 2 - 2Pv 3 

V\ 

\ Pv 4 - v 4 J 



w 2 - Pw 2 
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1 

+ 2 



+ 



(Pv 2 — Vi + v 2 \ 



2Pv 2 



/ ~2Pv 4 \ 


-2Pv 4 

Pv 4 + V 4 

ypw 2 + V 



W3 + PM/3 1 

* h - 

2 2 



U>4 + PlVi 1 

* h - 

2 2 



/^i - Pv 2 - v 2 \ 



-2P*; 2 



V 

/ 2Pt; 4 \ 




w 3 - Pw 3 



2Pvi 

Pv 4 — f 4 
\Pv 2 -v 2 J 



W4 — PW4 



where * is convolution. As required, T 2 (v, w) = T 2 (w, v) and T 2 (v, w) G W. 
For later reference (equation may continue on the next page): 







1 


'Pv 2 + v 2 — Pv 3 - v 3 ^ 












ssr 2 (v,w) = - 


vi - v 2 + Pv 2 
—Vi + v 2 - Pv 2 


* W\ 


V 1 


fvi + Pv 3 -v 3 \ ( \ 


1 
+ 2 


2v 2 + 2Pv 2 
—v 2 - Pv 2 


w 2 + Pw 2 1 

* h - 

2 2 


-Vi + Pv 2 - v 2 


u> 2 - Pw 2 

* 

2 


\ Pv 4 + V 4 J \ Pv A — Vi J 


(~Vi\ (-v 2 - Pv 2 \ 


1 

+ 2 






w 3 + Pw 3 1 

* h - 

2 2 






u> 3 - Pw 3 

* 

2 


\o ) \ ) 


(v A - PVi\ (v A + Pv A \ 


1 

+ 2 




PV4 + V A 


W4 + Pw 4 1 

* h - 

2 2 




Pv A — U4 


* W4 " PW4 (23) 


\Pv 2 + v 2 J \Pv 2 - V 2 J 



The nonlinear operator Q + : (0, 00) x Vcaugcd _ > W. Define 

n + (A,v) = SH A v + Ar 1 v+»r 2 (v,v) (24) 

with components denoted Q + = (Qf ,Q^ ^ti^t^t)- By construction, 

ToSeries(fi+(A, v)) = fi+ ction 2 2 (A, ToSeries(v)) (25) 

where ToSeries is applied component by component. Note that ^g ection 2 2 
redirects to (11), but that (11) must be interpreted as in Section 2.2. 
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Periodic 2D problem for Fourier-Chebyshev series (2DprobSeries). 

Find a solution (fi,co) € (0, oo) x Vcauged to fl + (fi,u) = 0. 

By construction, in particular (25): 

Claim. If (fi,u) solves 2DprobSeries, then the pair (/i, ToSeries(o;)) ; as 
a map on 3 section 2.2 for some £* > 1, solves 2Dpre and 2Dprob in Section 
2.2, except for (b) and (d), which must be checked separately. 



Remark. It is convenient to split Q + = into two parts. We have: 

SO + = and S fl fi + = <^ Q + = 

The direction ==>■ would be immediate from the definitions of S and S", if 
H" 1,re9 had a left inverse. Nevertheless, the following observations yield =^: 

• The first and third terms on the right hand side of (24) are in the image 

of E, and SS _1 ' reg S = 5. 

• The components 1,2,4,5 of i\v are all in the image of 1 — P, and 
5H" 1,teg (l — P) = (1 — P). Component 3 of i\v is annihilated by S. 



The (j system. Introduce the V — > V linear operators 

4» = diag^) 

Jjj, = diag(K M , J M ) = TauDerivative +/x jS* 
'0 



r w 

and the Y\ : V x V" — > V bilinear operator 

r j / \ l ( v 2 + Pv 2 ~ PV3 ~ V 3 \ 

r v.w = - „ } *w 1 

2V ' ; 2 V v i -v 2 + Pv 2 J 

1 f -vi + Pv 3 -v 3 \ ^ w 2 $ + Pw 2 $ 

2 \-v 1 + v 2 + 



■wi,)* 



1 (v 2 + Pll 2 + PV 3 + V 3 \ W 2 $ - Pw 2i 



2 V -vi + v 2 + 3Pv 2 J 2 
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where, v = (vi, V2,v 3 , u 4 ) and w = (w 1 ,w 2 $). 

Claim. If (fi,u) G (0, oo) x Vcaugcd is a solution to SQ + (fi,u) = 0, then 
u* = S tt r2 + (/i,w) satisfies the linear homogeneous identity 

J^J + fiE- l ' reg T\J + V\(u,J) = (26) 

See (16). Note that ToSeries(w tt ) = (^i, ^ 2 j) Section 2.3 by (25) and Sfi + = 0. 

4 Approximate solution yields true solution 

This section contains the analysis. We rearrange SQ + = 0, solve it with 
the contraction mapping principle, and then use the Jj system to show that 
Q + = 0. In other words, we solve 2DprobSeries. The calculations are 
modulo definitions (e.g. choice of explicit parameter values) and assump- 
tions (e.g. computer assisted results) that we postpone to Section 5. For 
convenience, things that we postpone to Section 5 are tagged by *. 



Cutoff operator. Fix integers m(s) > 2 and n(s) > 1 * 1 . Define s : V — > V 
and s : V -» V by s = diag(s, s,s,s) and 

J v mn if \m\ < 771(5) and \n\ < n(s) 
[ otherwise 



Reference. Fix a reference (/i rc f,u; rc f) G (0, 00) x Vcauged * 2 - Set 

(/!, U) = (/i rc f, CJ re f) + (/icorr, ^ corr ) (27a) 

This is a change of variables between the unknown (/i, oo) and the new un- 
known (/icorr, w corr ) G (— yU rc f , oo) x Vcaugcd- Both will be used simultane- 
ously, as convenient, with the understanding that they are always given 
in terms of one another by (27a). The condition fi co „ > —fi m { must be 
kept in mind. It is convenient to split the reference field itself into two 
parts, i.e. to fix (fi rclA , w rcfA ) G (0, oo) x V GaU gcd * 3 and to fix {fj, TeiB , Urem) e 
(-AWA, oo) x Vcaugcd * 4 , and to set 

(/i re f , U! re {) = (/U re fA, WrefAj + (/^refB, WrefBj (27b) 

Informally, refA is the main part, refB is for fine tuning. The calculations are 
organized in such a way that the most time consuming computer calculations 
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only involve refA. If the reader wanted to rigorously construct more accurate 
approximations to the Choptuik solution than those stated in this paper, he 
or she could use the same refA, but a better refB. The assumption /U re fA > 
implies that J^^ and J^ fA are defined. For convenience, suppose 

(l-sK cfA = * 5 (28) 

In particular, w rc fA has only finitely many nonzero components. 

Rearranging Sf2 + (/i, w) = 0: informal version. Set C = SQ + (fi,u) — 
Sf2 + (/i rcf ,u; rc f), and use the replacements (27a) and (27b) to exhaustion. Use 
the distributive law to multiply everything out. Clearly, C is quadratic in 
(/x corr , CJ corr ) without constant term. Split C = C\ + C 2 where C\ are all those 
terms in C that satisfy both of: 

• They are homogeneous of degree 1 in (/i CO rr, w corr ). 

• They do not involve (/x r efB,^refB)- 

All the remaining terms go into C^. At this point, write 

Ci = SQ + (fI,Uj) - Sfi + (/i r ef,Wrcf) ~ C 2 

By construction, C\ can be decomposed as C\ = 3^ IcfA uj corr + fi corr a + Au corr , 
where a is a multifield, and A is a linear operator. For high frequencies the 
term J MrcfA dominates. Therefore split A into the low-to-low frequency part 
sAs, and the rest A — sAs, and move the latter to the right: 

SQ + (fi, u) - Sfi + (/i rcf , Wref) - C 2 - (A- sAs)u co „ 

We obtain a rearrangement of SQ + = that will be useful to apply the 
contraction mapping principle, namely: 



Sn + (ji,u) = 



uj 



rr + UtetA V^ corra ~"~ ^ASUcquJ 

J MrcfA ( ~ Sfi + (// ro f, W re f) - C 2 - (A - sAs)u c 
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y corrJ 



J; ^corrj 



Rearranging SQ + (fi, u) = 0: rigorous version. Let C,Ci,C2,a G V and 
the linear operator A : V — > V be given by: 

C = Jfj, rcf W C orr + /^corrC^ +Sri)d; rc f + jl rc { SFiU; corr + 2S 5 ^(Wrcf, W corr ) 
+ ^corriji +Sri)d; corr + S^ r2(w C orr, ^corr) 
C\ = J^ refA UJ co „ + /i CO rr(4 +Sr!)u; rc fA + /ircfASFiOJcorr + 2S C r 2 (u; rc fA, UJ C 
C% = /i re fB A W corr + ^corrC^ + SFi)u; r cfB + ^rcffiSFiWcorr + 2S 5 ^(WrcfB, ^ c 

+ /^corrC^ +Sr 1 )u; corr + S^ r 2 (w corr , CJ corr ) 

a = (J+Sri)a; rc fA 

^4 = ^ rc fASr! + 2S S r 2 (a; r efA, • ) 

(29) 
It is now convenient to introduce the linear map (there are several equivalent 
ways to write this definition, because a = sa and J" 1 s = s J" 1 , s): 

J MrcfA MrefA / 

U : R © Vcauged -> V 

A © v ^ v + J- 1 ^ (Xa + sAsv) 

and the nonlinear map it : R © Vcauged — > V by 
*(A,v) = 

_J /^efA S ^ (^ref,Wref) 

~ J awa {/ i rcffi/<Iv + A(4+Sri)o; rc fB + /i rc f B Sriv + 2SHr 2 (a; r efB,v)} (30) 
-J-^jA^H-SrOvH-SHr^v)} 



MrcfA 

r 1 

MrefA 

r 1 

MrefA 

By construction 



S^ + (/i,w) = ^=^ U(// co rr © W corr ) = *(/i C orr,W corr ) (31) 

Remark. Domain and range of U decompose as 

(R© (VcaugcdH image s)) ©image(l-s) ->■ (images) ©(image(l-s)) (32) 

Note that image(l — s) C V Gauged by m($) > 2 and n(s) > 1. By con- 
struction, U is block diagonal with respect to this decomposition. Since the 
high-to-high frequency block image(l — s) — > image(l — s) is the identity, 
the operator U is invertible if and only if its low-to-low frequency block 
R © (Vcauged H images) — > images is invertible, a finite square matrix. 
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Weighted i 1 norms. Fix constants «i, k 2 > 1 * 6 and k* > * 7 . Set 

||z||c = I Re 2 1 + | Im^l 

||„,|| V~* /",, \\m\( ^. \\n\\\„, II 

IMIV - Z^m.neZ^lJ V^J Fmn||c 

= Em,„>o( 2 ~ ^mo)(2 - cU)K) m (>2) n ||w mra ||c 

ll v llv = IMIv + IMIv + ||^3||v + |K||v 
||A® v|| RffiV = K*|A| + ||v|| v 

with v = (wi, V2, vs, v±) G V. Note that || • ||c is a norm on C, viewed as a 
real vector space, and \z\ < \\z\\c < y/2\z\ and H-Zi^Hc < 1 1 -^l 1 1 <c 1 1 -^2 1 1 c ■ 

Note that || • ||v is not a norm on V, but only on {v E V | ||i> ||v < 00} C V. 
We do not introduce a special name for this subspace. 



Operator norm. For every linear operator O : X — > Y set 

||0||y<_jsf = SUp { ||Ox||y | X G X with \\x\\x = l} 

Then ||0'0||^ x < ||0'||^y ||0||y<-x. 

Convolution estimate. The estimate ||w * w||y < ||f ||y llvllv follows from 

\\V * W\\ V = E m ,nez( K l) H ( K 2) |n| || Em 1 +m 2 =m,n 1+ n 2 =n^m 1 n 1 W m2 „ 2 || c 

— 2-<m,neZV K l/ ^ 2 / 2^mi+m 2 =m,ni+n 2 =n ll^mim || C |Fm 2 n 2 II C 
^ Z^m 1 ,ni,Tia,na6ZV«l) 1 '(^J 11 2I ||Vm lie |Fm 2 n 2 lie 

— Z^mi,ni,m 2 ,n 2 GZV ft ' 1 / K^ 2 ) II y mim ||<D || ^Jm-xn-x ||C 

Here, H^i-^Hc < Ikillcll^Hc an d the assumption tti,«2 > 1 are used. 
This convolution estimate is sharp, in the sense that ||(w*)||y^v = IM|y. 
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Estimate for T 2 . 



|ssr 2 (v, -)llv<-v 

(Pv 2 + v 2 - Pv 3 - v 3 \ 
1 Vi — v 2 + Pv 2 

4 -vi + v 2 - Pv 2 



< max 



v 



(v\ + Pv 3 - v 3 \ 

2v 2 + 2Pv 2 

-v 2 - Pv 2 

\ Pv 4 + v A I 



V 



( 







—Vi + Pv 2 - v 2 
Vl 

\ Pv A -v 4 ) 
(-v 2 - Pv 2 \ 



\ 




1 


f-vi\ 




J 


V 


2 




I o J 



V 



V 









/ 



v 



(v A - Pv 4 \ 



PVi + v 4 

\Pv 2 + V 2 J 



V 



fv 4 + PvA 





PV4 — W 4 


\Pv 2 -v 2 ) 



(33) 



To prove this estimate, use (23) and 

||v||v = |K||v + \\\^2 + PvAv + |||^2 - Pv 2 \\ V 

+ |||W3 + Pv 3 \\v + |||^3 - Pv 3 \\v + ||K + Pv4\v + |IK - Pv4\v 

As a corollary, 



|SHr 2 (v, w)||v < 3||v||v||w||v 



(34) 



Estimates for the auxiliary operators U and V. 



U[fc,/]||v<-v< ( K a) fe sup mjri > ||/ mn || c <Cf(n 2 ) 

v[kj\\\ v ^ v <c f {i-{K 2 y k )- 1 



SU Pm,n>0 I/" 



with Cf — 1 if / m „ e R for all m,?7,, otherwise C/ = a/2. 

For V we assume, as always, that /c > and \f mn \ < 1 for all m,n. 



(35a) 
(35b) 



Proof: Let 



(2 - <S m0 )(2 - <5 n0 ) for all to, n > 0. Then 



U[fe,/]u||v = Em : n>0 C ™«( K l) m ( K 2)"||/mr 1 Wm,n+/c||c 

< ((K 2 )- k SU Pm „ ||/ mn || C ) Em,n>0 C ™,«+fc( K l) m ( K 2)™ +fc ||Wm,n+fc||c 



<((«2 



,-fe 



sup r; 



l/ mn || c )ikliv 
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Here c mn < c m! „ + fc for all m, n, k > has been used. Now 

II V[k, /]||v<-v = II E"=o U[fci, /l]||v<-v 
<Ei=ol|U[fcL,/ L ]||v^v 

< Ei=0 C fi,{K2)~ kL su Pm,« |(/i)mn| 

< c f (i - fa)-")- 1 

with (/ L ) m „ = {-l) L ]]%Io f m ,n+kt- Here C/ L < C f and |(/L) m „| < 1 have been used. 

Estimate for SI\. By (35), 

HSrjv^v < 4|| U M V 2il || v <-v < £! = 4(^)^(1 - («: 2 )" 2 )- 1 (36) 

To check this, note that the first component of STiV vanishes, for all v G V. 

Auxiliary estimates for U -1- ^. By (35), 

II J^ cfA ||v<-V < max{|| K^ fA ||v^ v , || J^ || v ^ v } 

< /C 2 d = f 2/^(1 - M" 1 )"^ + (/to)- 1 ) (37a) 

II J ^ fA ^llv^v</C 2 (37b) 

For (37), note that 

(i - (k,)- 2 )-^ + (^r 2 ) < a - (« 2 )- 1 )- 1 (i + (^r 1 ) (38) 

For (37b), note that J;V = ^ V[l, /](U[0, ^gg^ ]+U[l, W2 ff (w+1) ]) 
for some / with |/ mn | < 1. Similar for K M 1 / K. 



Suppose U is invertible * 8 . Fix an invertible linear operator 21 : V — > V 
with 21 = s2ls + (1 — s) * 9 . (The operator 21 is only introduced for an 
important operator norm estimate of the form HO1O2II < ||Oi2l~ 1 || ||2102||> 
to improve over the naive < ||Oi||||0 2 ||. See the left hand sides of (39d) and 
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< s 1 * 18 


(40a) 


< s 2 * 19 


(40b) 


< s 3 * 20 


(40c) 


< s 4 * 21 


(40d) 


<s 5 * 22 


(40e) 



(39e) below.) Fix constants < Ci, . . . , Cq < oo * 10 . Suppose 

||S~r 2 (u; rcfA , -)||v<-v< Cx * 1X (39a) 

Ikefsllv < C. 2 * 12 (39b) 

||Sfi + (/i rcf ,u; rcf )|| v <£3 * 13 (39c) 

||aJ^(A -*A>)||v«-v < A * 14 (39d) 

HU-^-inev^v^A * 15 (39e) 

HU-l^ev^v^A * 16 (39f) 

Then (this estimate is used only in Section 4.1) 

def 

||^||V«-V < ^3 = l^rciA^l + 2£i 

Fix constants < <Si, . . . , £5 < 00 * 17 . Suppose 

£ 6 /C 2 (l + /Ci)£ 2 /^* 

£ 6 /C 2 (/l rc fB + /^rcfB^l + 6£ 2 ) + £5/^4 

£ 6 /C 2 (l + K\)/k* 
3£ 6 /C 2 
Set iS 2 3 = max{5 2 , 1S3} and 1S45 = max{iS4, 25s}. 

Estimates for U -1 -^-. We have 

||U _1 *(A,v)|| ]RffiV < Si + S 2 n*\\\ +53||v||v + 54K*|A|||v||v + 5 5 (||v||v) 

< <Si +5 23 ||A© v|| RffiV + |<S 45 (||A©v|| ] a©v) 

(41a) 
by (30), (34), (36), (37), (39), (40). Furthermore, 

||U- 1 *(A 1 ,v 1 )-U" 1 *(A 2 ,v 2 )|| RffiV 

< \S 2 a + |5 45 (||Ai © v 1 || ]ReV + ||A 2 ffi v 2 ||]R.ev)J||Ai © vi - A 2 © v 2 || ]Rffi v 

(41b) 

Solution to Sf2 + (/i,w) = 0. Fix a constant 1Z > * 23 . Suppose 

S t + S 23 n + \S^TZ 2 < U * 24 (42a) 

S 23 + S i5 TZ < 1 * 25 (42b) 

Tl < K*/i rcf * 26 (42c) 
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Then Affiv4 U _1 *(A, v) maps B ->■ B with 



£ 



A©V G R©V G auged 



|A©v|| ReV <^} 



and is a contraction on the nonempty, complete metric space (B, || • ||k.©v) 
by (41). Therefore, it has a unique fixed point /i corr © cj coit . Equation (42c) 
implies that /i corr > — /i rc f, as required. The corresponding (/i,c<;) satisfies 
Sn + (//,w) = by (31). Note that 



||o;||v < oo 
This follows from (28) and (39b), and from 

||/^corr © ^corr||lR.©V ^ ^ 



(43) 
(44) 



Rearranging the (j system. Introduce A\ U", "A"" : V — > V linear maps 

A« = /^? 1 «r{ + r!(w Bftl -) 

u» = 1 + (J^J-Wttfi 

*" = -(JL fA )" 1 {(^efB + /ico 1 ,.)(4 tt + s- 1 ' rcs r») +r»K cfB + a; corr , •)} 

-(jLj- 1 ^-*^*) 

(45) 
Here, s = diag(s,s) : V — > V . This is an abuse of notation, but it will 
always be clear if s : V — > V or s : V — > V . Since (/i, u;) is now a solution 
to SQ + (fi, u) = 0, it follows from (26) and the construction (45) that 

UV = •V (46) 

Note that U" = sU's + (1 — s). Therefore, U" is invertible if and only if the 
finite square matrix sU% : images — > images is invertible. 



Norm on (a subspace of) V . Set ||w|| v « = ||iui||v + Il w 2ttllv f° r & U 
w = (wi, w 2 $) G V . This is a norm on {w G V^ | ||w|| v it < oo} C V". 



Estimate for Yl 



l^(v, 



| v tt^ v tt < max 



1 / -Vi + Pv 3 - v 3 

2 \-v 1 + v 2 + 3Pv 2 



v 2 + Pv 2 - Pv 3 - v 3 

V l - V 2 + Pv 2 

1 



v« 

v 2 + Pv 2 + Pv 3 + w 3 
— Vl + v 2 + 3Pv 2 



v : 



(47) 
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Auxiliary estimates for (U") l jcK We have 

I|r5,(v, -)|lv«4-v" - 3 ll v llv 
II H- 1,W8 rl|| v »^ v , < k\ = 2(K 2 y\\ - (k 2 )- 2 )- 1 

ikjLj-Iv**-* < *i = 2^a(i - (^)- l )- l (i + m- 1 ) 

WVLJ- 1 A* \\v^v« < Kl 

Suppose U" is invertible * 27 . Fix an invertible linear operator 2l" : V — > V 
with St» = <s2l«s H- (1 - s) * 28 . Fix < c\, 4, 4, 4, 4 < °° * 29 - Suppose 

||4(o; refA , OIU-v*^! * 3 ° ( 48a ) 

H^refB + ^corrllv < ^2 (48b) 

II^CJL.)" 1 ^ -*^«) IN*-* < 4 * 32 (48c) 

||(U»)- 1 (2t«)- 1 || v ^ vB <4 * 33 (48d) 

ll(U»)- 1 || v ^ v »<4 * 34 (48e) 



Then 



V tt^V» < ^3 = /VfA^i + 4 



Estimate for (UK) -1 *". 

||(U I, )- 1 * I, || V I. < _ V « < 1/irefB + /W|4*4(l + $) + ZC\lClC\ + 44 (49) 



Suppose the right hand side of (49) is < 1 * 35 . Then ||(U fl )- 1 * tl || v ^ v « < 1. 
We show that this implies cu" = 0, i.e. that (/x, u;) solves 2DprobSeries. 

The claim seems to follow from (46), which yields (1 — ||0|| V 3^v«)ll a; ''llv tl — with O = 
(U*) -1 -^ - ', and which seems to imply that ||w"|| v j must vanish. This reasoning is invalid, 
because ||w*|| v « is not known to be finite. We have to work a little more. 

For every v with max{l/«i, I/K2} < v < 1, set: 

IMI(v^) = E ro ,„>o( 2 ^ S m0 )(2 - S rM )(iyK 1 ) m (iyK 2 ) n \\v mn \\ c 

Define || • |j(v,^) an d || • ||(v".i/) accordingly. Note that || • ||(v,i) = II ' ||v- If f < 1, then 
(43) implies ||w'||(v« v \ < 00. Therefore, if we can show that ||0||(v*,i/)<-(v lt ,i') *"• ^ ^ or 
some v < 1, then w" vanishes. Since ||0|| V it^ v « < 1 by assumption, it suffices to prove 

limsup^n l|0||(v»,i/)<-(vV) ^ ll°llv»<-V« ( 50 ) 

36 



To check this, let e = (ei, e2j) G V be a unit vector, i.e. there are d{e) G {1, 2(1} and fc(e) G 
{0,1} and m(e),n(e) > such that (e d(e) ) m(e)j „ (e) = 1 if k(e) = and (e d(e )) m ( e)] „( e ) = i 
if k(e) = 1, and all other components of e are zero. Unit vectors are useful, because 

||Oe||(y|t v \ 

\\0\\(\;»^)^(\;»^) = sup .. .. — — (if the left hand side is < oo) 

e: unit vector II ^11 (V",f) 

This is a special property of the I 1 operator norm. Fix a unit vector e. For TV > 0, 
decompose Oe = (Oe) N+ + (Oe) N - with (Oe) N+ , (Oe) N ~ G V i and 

,,,,„», . n un \N+\ j((Oe)d)mn if m + n > m(e) + n(e) - N 

dG{l,2(t}, m,n>0: ((Oe) d + ) m „ = < 

I otherwise 

Suppose vq 1 v satisfy max{l/Ki, I/K2} < ^0 < v < 1- Then 

||Oe|| (vV) < \\{Oe) N+ \\ (ytM) + \\(Oe) N -\\ {vKu) 

< v m ^ +n ^- N \\{Oe) N+ \\ v t + {v/v ) m{e)+n{e) - N \\{Oe) N -\\ (v t :Uo) 

< v m(e)+n{e) - N \\Oe\\ vi + (^/^o) m(e)+Tl(e) ~ JV !|Oe|| (v ^ o) 

Since ||e|| (vM = u m ^+ n ^\\e\\ v , = {v/v Q ) m(e)+n{e) \\e\\ {vKvo) , it follows that 
l|Oe|| (vM N \\04 yL N l|Oe|| (vVo) 



ll e ll(V"^) ll e llv« l|e||(v»,r/ ) 

Since this holds for all e, and since ||0||/ v t v \^.(yt „) < 00 (see the remark below), 

l|0||(VV)<-(V«,r/) - ^"^H^llv^v" + { v / v oY N \\0\\(VKv )*-(vKv Q ) 

Fix i/ , set N(v) = (1 - v)- 1 / 2 and let w\l. Then v~ N M ->• 1 and (v / V ) _ " A ' M -» 0, and 
inequality (50) follows, because ||0||( V j „ \ w v « UQ \ < 00 (see the remark below). 

Remark: We show that ||0|Lyj v \<_r V » v \ < 00 if max{l/fti, I/K2} < ^ < 1. Separately, 

H('U 1, ) _1 ||(vit,i,) < -(vV) < °° ll^llcv.^MV,*) < °° 
The first by (U 8 )" 1 = s(U") _1 .s + (1 - s). The second by direct inspection of (45), using 

||WrefA.||(V,i>) < ||WrefA||v < °° and ll^rcfB + ^corr|| {V,u) < ||WrefB + W C orr||v < OO, See (43), 

and ll( J LfA) _1 ^ J H(vV)-KVM < °°> amon g others. 

4.1 Exterior estimates: (39d) and (48c) 
Fix cutoff operators t, q * 36 similar to s. Suppose 

2 < m(s) < m(q) < m(t) - m(s) + 1 < m(t) * 37 (51a) 

1 < n{s) < n(q) < n(i) - n(s) + 1 < n(i) * 38 (51b) 
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and 3 < n(q) * 39 . It follows that sq = qs = s and qt = tq = q and 
st = is = s. For convolution, q((t-) *(«■)) = and qSEr 2 (t-,s-) = 0, 
where a dot is any argument. In this section only, we abbreviate 

5 = 1 — 5, s = l— s, etc (the bar is not for complex conjugation) 



For every linear operator O : V — > V, the i 1 operator norm satisfies 

||0||v<-v = max{||Ot||v^v, ||Ot|| v <-v} 
Similar for V . Suppose the t parts in (39d), (48c) are under control: 

||*J;^(jl-*j4*)t||v«-v < U * 4 ° (52a) 

II^CJLJ- 1 ^ -*^)t|lv^v» < 4 * 41 (52b) 

Suppose that in addition, 

||2l||v«-v(£2.Mi + M 2 IC 3 ) + M 3 IC 3 < U * 42 (52c) 

Plv^v" (KlMi + M 2 K§ + M 3 K\ < C\ * 43 (52d) 

with 

Mx = 4/i rcfA ( K2 )-("«-"W +1 )(l - {k 2 )- 2 )- 1 
M 2 = 4/i r -j A ( K2 )-^-^ +1 )(l - M" 1 )" 1 
M 3 = 2(1 - (n 2 y i y l max {2 (m(q))" 1 , ^ A (n(q) - l)" 1 } (l + (k 2 )~ 1 ) 

In the rest of this section, we show that (52c), (52d) control the t parts. 
In other words, we show that (52) implies (39d), (48c). 

For (52c), set O = QIJ~ (A — sAs), a local notation. Since st = 0, one has 

It follows that 

||Ot|| < ||2t|| || J7* || ||qAt|| + ||2t|| \\sJ-\ A q|| ||A|| + || J;\ A q|| ||A|| 

II II — II 1 1 1 1 /l r cfA nil' II II 1 1 1 1 MrcfA ' 1 1 1 1 II H /J rc f A ' 1 1 1 1 H 
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with || • || = || • ||v«-V> a local notation. Recall (28). We have 

||qAt|| = /x rcfA || qSrxtH < /i rcf A4||qU li iV 2i it||v^v < Mi 

II s J ^efA All ^ max ill s ^L *\\v<-v, HsK-^qllv^v} < M 2 
II J- 1 ^ q|| < max{|| J" 1 ^ q|| v <_ v , || K" 1 ^ q|| v ^v} < M 3 

Therefore, (52c) implies ||Ot|| < £4, as required. 



In (A), wc estimate H-SJ^^ q||v<-v- Similar for ||sK^ fA q||v<-v and ||q Ui,i V2,i t||v<-v- 
In (B) we estimate || K~ rsfA q||v-<— v- Similar for || J^ refA q||v<-v- 

(A) There are numbers Q mnn ' e € such that (s J^ fA qu)mn = E n '>o Qmnn'V m n' for all 

V G V, and ||Q mnn /||<rj < 4 MrefA' Als0 

Qmnn' 7^ ==> ( (m < m(s)) and (n < n(s)) J and ( (m > m(q)) or (n' > n(q)) ) 
Since m(s) < m(q), this simplifies to 

Qmnn' 7^ => (m < m(s)) and (n < n(s)) and (n' > n(q)) 
Therefore, with c mn = (2 - 5 m0 )(2 - S n0 ), 

ll S 3 »LaW\\v = J2 m ,n>0 C ^n(Kl) m (K 2 ) n \\ E„'>0 Qmnn'Vmn' || c 

= Lo<m<m(»),0<n<»(i) Cmn(^l) ( K 2J || En'>ra(q) Qmnn'Vmn' \\ c 

— ^rcfA Lo<ro<m(s),0<n<n(») Cm»(«l) (ts) En'>n(q) IKrm'llc 

^ 4 A«refA2_/0<m<m(s),n'>n(q) C ™«'( K l) m ( K 2)" || «mn' || C Eo<n<n(«i) ( K s) ™ 

<4/X r -J A (E m ,„'>0^n'(«l) m («2)"'||t;w|| C )(Eo<n<n( S )(«2) ( "-^ )) ) 

< v-fAiiHiv^r^- 1 -"^^! - m- 1 )- 1 

Here c mn < c mn < has been used in one step. 

(B) Let r be the cutoff with m(t) = m(q) and n(r) = ?i(q) — 2, a local notation. Then 
K^ fA q = V[2, /] L- 1 ^ (1 + LVOq = V[2, /] L" 1 ^ f (1 + U 2 ,_i)q for some /. Therefore, 

II K »L qllv^v < II V[2, /]|| v <-v|| L-^ [A t|| v ^v(l + II U 2 ,_ x || v ^ v ) 

Use || V[2, /]|| v <-v < \/2(l - M -2 )" 1 and 1 + || U 2 ,_i || v <-v < 1 + M" 2 and 

||lL-t fA t||v^v<V^ sup ^- i X -^<V2max f 



below |im/2 + /UrefA(n + 1)| I m(q) /irefA(n(q) - 1) 

where the sup is over all m, n > with to > m(q) or n > n(q) — 2. Now use (38). 
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For (52d), set O i = ^(Jj^j- 1 ^ - sAh). Then 

\\o r t\\ < W\\ ikjL a ) _1 II IM*tH + Ml \WLJ- X *\\ UH 

with || • || = || • Hv^v"- R eca ll (28). One has 

||qA*t|| =/i rcfA ||qH- 1 ^r tt 1 t|| </i re fA2||qU 1:1 V 2il t|| v ^v<Mi 

\wLJ-^\\ < m 2 

ll(JLJ^H<^3 

Therefore, (52d) implies ||0"t|| < £ 3 , as required. 

5 Parameter values and Computer assisted results 

This section takes care of every individual * tag in Section 4. Set 

Kl = 65/64 k 2 = 5/4 k* = 1 

m(s) = 250 m(t) = 1600 m(q) = 1350 

n(s) = 750 n(t) = 4800 n(q) = 4050 

We have * 1§/ * 6 " * 7,/ * 36 " * 37 " * 38 " * 39 ". The text files RefA.dat and 
RefAplusB.dat contain (/i rc fA? ^rofA) and (/x re fA + /^refB, ^rcfA + Wreffl); the file 
format is self-explanatory from the discussion in Section 6 * 2 * * 3 * * 4 ^ * 5 1/ . 
One finds, either by looking at the first few lines of these files, or with the 
command choptuik LoadAplusB: 

/i rcfA = 722873400 ■ 2^ 32 
/WA + /ircfB = 786320438575298 [. . . 180 final digits omitted . . .] • 2" 650 * 

|/ircfB|<43-2- 40< ' 

Most readers should ignore the "* tags. (They tag the things that must be 
revisited if one repeats the construction with a different refB, but the same 
refA, to get a better approximation to Choptuik's solution.) Set 



£1 = 10.5 


£ 2 = 2" 25 * 


£ 3 = 2~ 294 " 


£ 4 = 0.625 


£ 5 = 1.004 


£ 6 = 258 


Si = 2~ 279 * 


S 2 = 0.01* 


S 3 = 0.64* 


S 4 = 2 19 


S 5 = 2 17 




4=8 


4 = 2~ 25 * 




4 = 0.625 


4 = 1.004 


4 = 258 



40 



*!()• *17«/ *29t/ 



(»v) 



Set 21 = diag(03, 53, 23, 03) and 21* = diag(03, 03) where 

256t> mn if m < 50 and n < 150 

16u mn else, if m < 100 and n < 300 

4w mn else, if m < 150 and n < 450 

2w mn else, if m < 200 and n < 600 

1.5625w mn else, if m < 250 and n < 750 

^?nn else 



in particular, 03 = s03s + (1 - s) * 9,/ * 28t/ . It follows that 

/d = f /C 2 < 107 /C 3 < 22.5 

<S 23 = 0.64" 545 = 2 19 

K\ = f K\< 107 /C* < 8.75 

M x < 2" 240 M 2 < 2~ 1055 M 3 < 0.027 

PHv^v < 256 Plv^v" < 256 

j^ QW *18 • *19 • *20 • *21 • *22 • 



(39a): Use (33) with v = w rcfA . Use choptuik DedicatedGAMMA2 . 

(39b): Use choptuik LoadAplusB. ** 

(39c): Use choptuik FailureToBeSolAplusB. "" One output line is 

(11 norm of SBold(0megaPlus(ref AplusB))) <= 15312135*2" (-320) 

(This command requires ~ 16 GB RAM.) This takes care of * 1:L * 12 * 13 1/ . 

We note that the command choptuik FailureToBeSolAplusB also yields 

(11 norm of OmegaPlus(ref AplusB) ) <= 24707998*2" (-320) 

which includes the constraint equations. We don't need this inequality. 

For (39d), note that (52c) holds. To check (52a) use the command 

choptuik EXTERIOR.ESTIMATE of fm_of fn_numm_numn (53) 

where of fm, of fn, numm, numn have to be replaced by actual numbers that 
satisfy of f m, of f n > and numm, numn > and offm + numm < m(t) and 
of fn + numn < n(t). For example, an actual command would be 
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choptuik EXTERIOR.ESTIMATE 212_748_2_4 

The command (53) generates a file with a number of lines, one line for each 
unit vector in V whose Fourier-Chebyshev indices (m, n) satisfy of f m < m < 
of f m + numm and of f n < n < of f n + numn. Each line begins with d_k_m_n, 
for the unit vector e = (e±, e^ e%, e^) G V with (e d+ i) mn = 1 if k = 0, and 
(ed+i)mn = i if k = 1, all other entries equal to zero. For example, the line 

1_0_212_750 550*2" (-10) TRUE (54) 

for the unit vector (e 2 ) 212,750 — 1 asserts that 

||St J~^ A (A - sAs)e|| v /||e||v < 550 ■ 2" 10 (55) 

The TRUE indicates that this number is < 0.625, as required by (52a). 

The file EXTERI0R_ESTIMATE_0_0_1600_4800 .dat was generated on an HPC 
cluster, and contains the result of applying (53) to all 26868000 unit vectors in 
image t C V. Many tools can be used to analyze the file. For example, every 
line of the file has the same structure as (54), because the unix command 

grep -v -c , ~[0123]_[01]_[0-9]*_[0-9]* [0-9] *\*2\~ (-10) TRUE $' 

EXTERIOR ESTIMATE 1600 4800.dat 



yields 0. Without the -v option, it yields 26868000. 



*14i/ *40«/ *42«/ 



For (39e) and (39f), use the command 

choptuik INV of fm_offn_numm_numn (56) 

with of fm + numm < m(s) and of fn + numn < n(s) (more details below). It 
shows that there exists a linear map V : V — > K, © Vcaugcd, block diagonal 
with respect to the decomposition (cf. (32)) 

( images) © ( image(l - s)) -)■ (R © (Vcaugcd l~l images)) © image(l - s) 

such that the high-to-high frequency block image(l — s) — > image(l — s) is 
the identity, and such that 

||1 - UVHv^v < 2- 16 llVSri^ev^v < 1 (57) 

(The operator 21 was constructed so that the second inequality would hold.) 
It follows that U is invertible, and given in terms of the Neumann series 

1 1 

IT 1 = ¥ ; r or IT 1 = V2T 1 ; -21 (58) 

l-(l-UV) 1-(1-21U¥21~ 1 ) v ; 
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The second Neumann series converges by 

||i - auvsr 1 !! = p(i - uv)^- 1 !! 

< ||2l|| ||1 - UV|| ||2t _1 || < 256 ■ 2~ 16 • 1 = 2~ 8 < 1 
The estimates (39e), (39f) follow from 

HU-^IUv^-v < (1 - 2- 8 )- 1 llU-lmev^-v < (1 - 2~ 8 )- 1 2 8 

The file INV_0_0_250_750 . dat was generated on an HPC cluster. It contains 
the result of applying (56) to all 654375 unit vectors e G images C V. For 
example, it contains a line for the unit vector with (64)1,7 = i, 

3_1_1_7 580*2" (-32) TRUE 1156600*2" (-13) TRUE 
This line asserts the existence of a vector / e R © (Vcauged H images) with 
||e - U/|| v /||e||v < 580 ■ 2" 32 ||/|Uv/l|e||v < H56600 • 2~ 13 (59) 

The first TRUE indicates that 580 • 2" 32 < 2" 16 , the second indicates that 
1156600 ■ 2- 13 < ||2te|| v /||e|| v , where ||2le|| v /||e||v = 256 for this particular 
e. Partially define V by Ve = /. Analogous statements for all unit vectors 
e define V completely, and yield (57). * 8l/ * 15l/ * 16|/ 

Set 

U = 2- 277 < 10- 83 * 1 (60a) 

*23»/ *24i/ *25^ *26./ rp^ resu iti n g fi xec [ p i n t (/i, w) satisfies 

|/i-(/ircfA + /ircf B )| <2- 277<< || LO - (w rcfA + W rcfB ) || V < ^^ * (60b) 



We proceed analogously with the Jj system. 

(48a): Use (47) with v = w rcfA . Use choptuik DedicatedGAMMA2. 
(48b): Use the triangle inequality, choptuik LoadAplusB and (60). ** 
(48c): Note that (52d) holds, and to check (52b) use the command 

choptuik SHARP_EXTERIOR_ESTIMATE of fm_of fn_numm_numn 
or see SHARP_EXTERI0R_ESTIMATE_0_0_1600_4800.dat 

*30t/ *31t/ *32t/ *41i/ *43t/ 
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For (48d) and (48e), use 

choptuik SHARP_INV of fm_offn_numm_numn 
or see SHARP_INV_0_0_250_750.dat 

There exists a linear map V" : V tt — > V tt with V^ = sV^s + (1 — s) and 

Hi - u»v«|| v ^ v » < 2- 16 \\vK^)-%^ v , < 1 

*27^ *33^ *34i/ ^ Qw *35t/ ^y g &re ^ one ^^ & l| ^e A ^ g ec tJ on ^ anc i 

we now have a solution to 2DprobSeries. 

To see that we also have a solution to 2Dprob, i.e. a solution to the Einstein 
equations with massless scalar field, we still have to check (b) and (d) in 
2Dpre with £» = \{k 2 + l/« 2 ) = 41/40. Let 3 be the cylinder (12). Then 

sup( T] £) G 3 | ToSeries(w)| < \\v ||v 

follows from \z\ < \\z\\c, and \e ±id \ < k 2 for all 9 G C with cos^ G (— £*,£*)■ 

Since ||/WB © W rc fB + /icorr © Wcorr||R©V < 2 ~ 24 *j We naVe 

I At - AWa| + St=i su P(r,c)e3 I ToSeries(wj - w RcfA ,j)| < 2~ 24 ** 
Now, (b) and (d) follow immediately by analyzing the data in RefA.dat. 

For (2) and (3), note that \ji — (aWa + AAreffi)! < 10~ 83 "* and 
K - ( 27r(/i rcfA + /ircfB) - / dr (w RcfA ,3 + WrbB^Ct, 0) J 

< 27r||/X cor r © W CO rr||lR.eV < 10~ 82 " (61) 



See (15b). Now, the data in RefAplusB.dat, or more conveniently the com- 
mand choptuik LoadAplusB, yield inequalities (2) and (3). 

6 Data structures 

To enable the interested reader to translate between this paper and the 
accompanying C source code in the directory sourcecode (arXiv ancillary 
files), we introduce the basic data structures that we use, and explain what 
data they represent. 

The type Sector contains four non- negative integers, 
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typedef struct 
{ 

unsigned long off_m; 

unsigned long off_n; 

unsigned long num_m; 

unsigned long num_n; 
} Sector; 

and represents a rectangular region of Fourier- Chebyshev indices (jn,n): 

of f _m < m < of f _m + num_m 
of f _n < n < of f _n + num_n 

The type GI contains two arbitrary precision integers (mpz_t from gmp.h): 

typedef struct 
{ 

mpz_t Re; 

mpz_t Im; 
} GI; 

and represents the Gaussian integer Re + i Im e 7L + i7L. The type 

typedef struct 
{ 

long Two Exp; 

Sector sec; 

GI** data; 
} Field; 

represents the field v = (v mn ) G V given by 

^sec.off_m+u, sec.off_n+v = 2 data[u] [v] 

for all u and v with: < u < sec .num_m and < v < sec .num_n 

and v mn = for all m,n > not in sec. The offsets off_m, off_n allow 
us to efficiently store and calculate with fields v that are 'strongly localized' 
in frequency space. For example, for tasks (2a), (2b), (3) in Section 1.3, we 
often use fields with of f _m ^> num_m or of f _n ^> num_n. The type 

typedef struct 
{ 

long TwoExp; 

mpz_t coeff; 
} DyadicQ; 
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represents the dyadic rational number coeff ■ 2 TwoExp . The type 

typedef struct 
{ 

unsigned long num_comp; 

Field* comp; 
} MultiField; 

represents the multifield (comp[0] , . . . , comp [num_ comp- 1] ) G y num - comp , and 

typedef struct 
{ 

DyadicQ mu; 

MultiField mf ; 
} mu_MultiField; 

represents (mu, mf ) 6R® V mf ■ num - com P. One can now also read the text files 
RefA.dat and RefAplusB.dat, each contains a mu_MultiField. 
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